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Abstract 

We derive and calculate unitarity bounds on the slope and curva- 
ture of the strangeness-changing scalar form factor at zero momentum 
transfer, using low-energy constraints and Watson final state interac- 
tion theorem. The results indicate that the curvature is important 
and should not be neglected in the representation of experimental 
data. The bounds can be converted also into an allowed region for 
the constants C[ 2 and C34 of Chiral Perturbation Theory. Our results 
are consistent with, but weaker than the predictions made by Jamin, 
Oiler and Pich in a coupled channel dispersion approach based on chi- 
ral resonance model. We comment on the differences between the two 
dispersive methods and argue that the unitarity bounds are useful as 
an independent check involving different sources of information. 
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1 Introduction 



The scalar form factor relevant for the decay K — > iriu is proportional to the 
matrix element (K\su\n). It is related to the vector form factors f±(t) by 
the expression 

/o(*) = M*) + M / M| /-(*), (i) 

in the standard normalization where /o(0) coincides with /+(0), a quantity 
that is of central importance for the determination of the CKM matrix ele- 
ment V us . Writing the Taylor expansion of the form factor near t = 

fo(t) =U(0){l + X t + ct 2 + ...} , (2) 

where A is the slope and c the curvature, it is customary to relate the slope 
with the scalar radius of the pion, through A = (r 2 )f n M 2 /6. This expan- 
sion is usually considered in the analysis of experimental data on decays 
of both charged and neutral kaons. The experimental situation on the Ki% 
decays improved very much recently. For the neutral kaons, the new result 
from KTeV collaboration is (r 2 )^ L7V = 0.165 ± 0.016 fm 2 [1. This value now 
dominates the statistics and is consistent with the 1974 high statistics exper- 
iment j2] • For the charged kaons, the data collected at the ISTRA+ detector 
give the radius (r 2 )f ±7r = 0.235 ± 0.014 ± 0.007 fm 2 , which dominates the 
world average. 

On the theoretical side there is also a considerable progress in the study 
of the scalar strangeness changing form factor fo{t). Low energy theorems 
impose strong constraints on the values of this function at some particular 
points. The value at the origin /o(0) ~ 1 is given by Ademollo-Gatto rela- 
tion 0], while the theorem of Callan and Treiman, refined by Dashen and 
Weinstein states that in the limit of vanishing quark masses the value of 
the form factor at the special point t = M\ — M 2 is 

f (A K7r )=F K /F n + A CT , A Kn = M 2 K -M 2 , (3) 

where the 0(rh) correction Act was calculated to one loop in Chiral Pertur- 
bation Theory (;\RT) jH] and found to be tiny: Act w —3 ■ 10~ 3 . Also, in 
Ref. a first prediction to one loop for the radius was obtained, (r 2 )^ n = 
0.20 ± 0.05 fm 2 . It is shown that the corrections of 0(rh) do not contain a 
chiral logarithm of the type M 2 log M 2 and therefore they are very small 6 a . 
Recently, the form factors were calculated to two loops of %PT 0IB]. 
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Dispersion theory provides another useful tool for making predictions on 
the form factors that govern the Ki 3 decay. These functions are analytic 
in the t-plane cut from the unitarity threshold (Mk + M^) 2 to infinity. A 
first approach to derive bounds on the values of the form factors and their 
derivatives in the physical region was proposed in the years 70. It is based on 
analyticity, and exploited the Kallen-Lehmann representation for a suitable 
two-point function if)(q 2 ), in combination with unitarity and low-energy the- 
orems. However, these first studies [ID]- [Hi were based on an unsubtracted 
dispersion relation for ip(q 2 ), which turned out to be incorrect: this function, 
calculated in perturbative QCD [T5] . satisfies a dispersion relation which re- 
quires two subtractions. The correct dispersion relation was applied for the 
first time in Ref. ^H]- Instead of using the low energy theorems on ip(0) as 
input PO]- [131 the new approach involves the calculation of the derivative 
ip"(Q 2 ) obtained by perturbative QCD in the deep euclidian region. The 
constraining power of this condition is weaker and therefore the strength of 
the bounds is reduced. 

An alternative dispersive approach based on coupled channel unitarity 
equations for the Kir, Krj and Krj form factors was considered recently in 
|17j . It is based on unitarized chiral perturbation plus a K- matrix fit of the S- 
wave scattering [TSj. In a recent paper [T^], the authors combine the values of 
the slope and curvature of the scalar form factor at zero momentum transfer, 
calculated in [T7j, with the two-loop calculation of |7J |S], obtaining 

thereby predictions for the chiral constants C[ 2 and C^. 

In the present work we revisit and extend the technique of unitarity 
bounds with the aim of deriving constraints on the slope and the curva- 
ture appearing in the Taylor expansion (J2J). The study is motivated by the 
recent interest in the determination of these parameters, both experimentally 
and in %PT. On the other hand, since the last work on unitarity bounds [TB] . 
the knowledge of the quantities entering as input improved considerably. For 
instance, the correlator ip(q 2 ) is now calculated in perturbative QCD to order 
°P S EB) E3, an d recent lattice calculations of the light and strange quark 
masses j2S], [231 an d the ratio F K /F n j2H| are available. We notice also that 
the implementation of Watson final state interaction theorem |2E! was never 
done with the correct dispersion relation for ip(q 2 ) (in Refs. J3]; where 
this problem was considered, the old unsubtracted dispersion relation was 
used). In the present work we implement simultaneously both low energy 
constraints and the Watson theorem, exploting in an optimal way the for- 
malism, which can be compared with the alternative study made in Ref. 
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In Section 2, we give a short review of the mathematical technique we will 
use, and show how to incorporate additional constraints provided by x?T and 
Watson final state interaction theorem. In Section 3, we discuss our results 
and compare them with the recent predictions on the slope and curvature 
of the scalar Ktt form factor presented in jlZ]-^!- As in these references 
we work in the isopsin limit (the isospin breaking corrections calculated in 
[2*Uj are important and should be included in a more precise analysis). We 
end the paper with some comments on the usefulness of the method and its 
possible generalizations. 



2 Unitarity bounds 

The starting point in the derivation of the bounds is the two-point function 
^j(q 2 ) = i f d 4 xe i<? ' x (0|T(9 M ^(x)9 M ^(0) t )|0) , (4) 



where = s'y^u is the strangeness changing vector current that governs 
Kiz decay. The function ip(q 2 ) was calculated for euclidian arguments Q 2 = 
— q 2 > up to corrections of order [T^], [21], |22j It grows asymptotically 
as Q 2 and satisfies a twice-subtracted Kallen-Lehmann representation. One 
can use either the dispersion relation for the second derivative ip {Q 2 ) 



oo 



t+ 

where t + = (Mr + M^) 2 is the first unitarity threshold, or the relation 
satisfied by the derivative of the ratio i[>(Q 2 )/Q 2 : 

{ ) + Q* ~ ttJ M t(t + Q*)f W 

The advantage of the last relation is that it incorporates the value ip(0), 
which satisfies the inequality 

<K0) < (M K F K - M^) 2 , (7) 

derived by Mathur and Okubo (we use the normalisation F n = 92.4 
MeV). In our analysis we shall use the relations © and (jBJ) successively. 
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Unitarity bounds ^U] are obtained by combining the above dispersion 
relations with the inequality 



3(Mt — M 2 ) 2 \ (t — t+)(t — t_) 
Im^)>^4- ^ -\fo(t)\ 2 , t ± = (M K ±M 7T ) 2 , 



64vr 



t 



(8) 

given by unitarity and the positivity of the spectral function. Then Eq. (|5|) 
gives [12] 



7T 



- dtw(t)\f (t)\ 2 <ip'\Q' 



where 



w(t) 



3(M| - Ml) 2 yj(t-t+)(t-t-) 



32vr 



Similarly Eq.© leads to 



oo 

-/dt^(t)|/o(t)| 2 < 

7T J 



t{t + Q 2 y 



>(Q 2 )V , V(0) 



where 



u)(t) = 



v Q 2 

t + g 2 

2t 



Q 4 



(9) 



(10) 



(11) 



(12) 



The r.h.s. of (jHJ) is given by perturbative QCD [T5].|22j 



3 (m s (Q 2 )-m u (Q 2 )) 2 

^7T 2 Q 2 



x 



11 « s (g 2 ) /Wq 2 

1 + + 141? _^v^ 

3 7T V 7T 



+ ... + — +0 — 



Q 



Q 1 



(,13) 



while the corresponding expansion for the r.h.s. of (jlljl reads 
/>(Q 2 )Y _ 3 K(Q 2 )-m„(Q 2 )) 2 



1 + 11 *m + 45.84 

3 7T V 7T 



Q 2 

a s (Q 2 Y 



m 



+ -.. + — +0 — 



Q 



(.14) 



Using standard techniques [ID] -[Hi], [2*7] . the conditions Q or (fTT|) can be 
converted into bounds on the values of fo(t) for points inside the analyticity 
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domain, in particular, on the slope and the curvature at the origin. We 
consider in detail the condition (jUJ), indicating at the end of the section the 
modifications resulting for the condition (fTTj) . 

The problem is brought into a canonical form by making the conformal 
mapping 



z(t) = V^-VgZ I, (is) 



which maps the t-plane cut for t > t + onto the disk \z\ < 1, such that 
z(0) = 0, z{t + ) = 1, z(oo) = —1, and the upper edge of the cut t > t + 
becomes the upper semicircle z = exp(i9), < 9 < it. The inverse of (I15j) 
and its derivative are 

t(z) = At + z/(l + z) 2 , dt/dz = 4t+(l- z)/(l + z) 3 , (16) 

while for points on the boundary, z = e ld , Eqs. (JT3J) and (fTtj)) write 

9 = 2 arctan ^*/*+ - 1, t(e ie ) = t+f cos 2 -6 . (17) 
In the new variable the relation (jSJ) becomes 

^ jd9p(9)\f(e W )\ 2 <^"(Q 2 ), (18) 
where /(z) = f (t(z)) and 

p(0) = «;(f( C *)) 



dt(e' 



■i0\ 



d9 



(19) 



We define the outer function C(z), analytic and without zeros inside \z\ < 1, 
such that on the boundary of the unit disk: 

\C(e»)\ =y/M- (20) 
A straightforward calculation ^3], [27] gives 



C( g ) -\r MK ~ M * (1 z)(l I z)^ V 1 ~ - + ^ + Z) (21) 
C[Z) - 4V2vr M^ + M/ 1 + (1-^ + ^(1 + ^)3' (21) 
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where (5 = yl — £_/£+ and /3q = J 1 + Q 2 /t + . Then the inequality (fTHj) can 
be written in the form 

^ /d%(e^)| 2 <V/'(Q 2 ), (22) 

Z7T J 

— 7T 

where the function g(z) defined by 

g(z) = C(z)f(z) (23) 
is analytic inside the disk \z\ < 1. Expanding this function as 



g{z) = £ g n z\ (24) 

n=0 

we write the L 2 norm condition (|22|) in terms of the Taylor coefficients g n as 

oo 

E^n<^"(Q 2 ). (25) 
n=0 

We mention that the reality condition fo(t*) = fo(t) means that the coeffi- 
cients g n are real. 

From Eq. (J23|) it follows that each Taylor coefficient g n is a linear ex- 
pression of the first n derivatives of the form factor fo(t) with respect to t at 
t = (we recall that t = is transformed by (|15|) into the origin z = 0), in 
particular g Q , g\ and gi are expressed in terms of the parameters /o(0), (r 2 )^ n 
and c. The coefficients depend on the conformal mapping z(t) given in ()15jl 
and the outer function C(z) given in (|21|) . where Q 2 enters as a parameter. 
For instance, for Q 2 = 4GeV 2 , a straightforward calculation gives: 

9o = /o(0) 0.00164, 

= /o(0) [-0.00189 + 0.0113 (r 2 )f^], (26) 
^ 2 = / (o) [-0.000259 - 0.03579(r 2 )f ^ + 0.00428 c]. 

In order to find the maximum domain allowed for the parameters /o(0), 
(r 2 )f v and c, we fix the first three coefficients go, g\ and #2, and take the 
minimum of the left hand side of (|25J) with respect to the remaining param- 
eters g n , n > 3, which are arbitrary. Clearly, the minimum is reached for 
g n = 0, n > 3, so the allowed domain becomes 

g 2 o + gl + gl<^"(Q 2 ). (27) 
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Using Eq. ([26)1. this inequality defines a convex domain for the parameters 
/o(0), (t" 2 )^ 71 and c. For a fixed /o(0), this domain is the interior of an ellipse 
in the plane ((r 2 )f ^ c). 

The strength of the bounds, expressed by the size of the allowed domain, 
is increased by imposing additional contraints. 

(1) The Callan-Treiman condition 

Let us consider first the Callan-Treiman condition (j3J). We note by zq the 

image of the point t = A Kn in the z-plane: z = V Mk / m ^^_ ^ From Eqs. 

and (}2~2j) we have the condition 

oo 

EM"-^o)/o(A fo )- (28) 

n=0 

As explained in Ref. [27], in order to obtain the allowed domain, we have 
to evaluate the minimum of the left hand side of (|23j) with respect to the 
coefficients g n , for n > 3 subject to the condition (j2*%|) . We impose this 
condition by the Lagrange multipliers method. The Lagrangian is 

oo / oo \ 

£ = E 9 2 n + I* E - C(z )/o(A iC7r ) (29) 

n=0 \n=0 / 

where /i is the Lagrange multiplier. The minimizing condition 

dC 

~ 0, n>3, (30) 



has the solution g n = —fiz^/2 for n > 3. Inserting this solution in the 
condition (}2*H|) . we find the Lagrange multiplier 



EHJi = 2 (g + 5l z + <?2*o - C(z )f (A K7T )) . (31) 



2i 







By inserting the solution of Eqs. (j3*Uj) with \i from (}3*9*|) into the left hand 
side of Eq. (|23jl . we obtain the inequality 

0o+0?+0i + ^— U + ^i^o + ^o-C(^o)/o(A A ^)) 2 <^"(Q 2 ), (32) 

which represents again, for each /o(0), the interior of an ellipse in the plane 

((r 2 )^, c). 
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(2) The Callan-Treiman condition and the Watson theorem 
We now include Watson theorem j2H]> which states that 

Arg[/ (t + ie)) = 5l /2 (t) , (M K + M w ) 2 <t<t in , (33) 

1 /2 

where 5 (t) is the 1 = 1/2 phase shift of the S wave elastic scattering 
and ti n — {Mk + 3M„-) 2 is the inelastic threshold. The incorporation of 
Watson theorem in the unitarity bounds for the Ki 3 form factors was done 
in an approximate way in Ref. and rigorously in Ref. [Tl]. However, 
as mentioned above, in these works it was assumed that ip(q 2 ) satisfies an 
unsubtracted dispersion relation, which was not confirmed in perturbative 
QCD. Moreover, the authors derived bounds on the magnitude |/o(£)| in 
the physical region, and did not consider the slope and the curvature in 
which we are now interested. In the present paper, we apply a mathematical 
method based on Lagrange multipliers, which is equivalent to the technique 
adopted in [13] but is conceptually simpler. The details are explained in Ref. 
[2~Tj . where the allowed domain for the coefficients g n with the constraint 
(|33|) was obtained. Below, we consider the complete problem by imposing 
simultaneously the conditions (jSJ) and (J33J). 

Let us note by 9 in = 2 arctan Jtin/t+ — 1 the image of the point t; m on 
the boundary of the unit disk \z\ < 1 in the z-plane. We define the Omnes 
function [Hj, £8! 

o{z) = ^jde T ^L g , (34) 

— TT 

where 8(9) = 5 1/2 (t(e^)) for < 9 < 9 in . When 9 > 9 m the function 5(9) is 
arbitrary, however these values will not enter the result. We recall again the 
reality condition fo(t*) = fo(t), which means that 5(9) is an odd function 
of 9: 5 (—6) = —5(9). As in Ref. [Hj, we assume that 5(9) is a Lipschitz 
continuous function of order larger than 1/2 on [0, 9 in }. By construction, 0(z) 
is analytic and without zeros inside \z\ < 1, and ArgO(e ?61 ) = 5(9). Therefore, 
in the ratio f(z)/0(z) the phases of the numerator and the denominator 
compensate, leading to a real quantity along the elastic region. Expressed in 
terms of the function g(z) defined in ()23)1 and its Taylor coefficients g n , this 
condition reads 

= 0, -9 in <9<9 m , (35) 



lim Im 

r-»l 



g{re 



i0 , 



W(re l 



lim^# 
i — >i — * 



Im 



n=0 



W(re 



> 8') 
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where we introduced the function 

W(z) =C{z)0{z). 



(36) 



We must find the minimum of the left hand side of (j2*5Jl with the constraints 
and (|3*5j). We consider the Lagrangian [27] 



£ = Y,9n + ME 9nzl - C(z )fo(A K7r ) 



n=0 



\n=0 



+ ]imJ2g n I ^\{0>)\W{e if>, )\lvy 



W(re i9 '] 



(37) 



where X(6) is a generalized Lagrange multiplier, which must be an odd func- 
tion: A(-0) = -A(0). As discussed in Ref. [27], the factor |t^(e ie )| was 
inserted in the integrand in order to have a function X(8) with integrable 
squared modulus (i.e. of class L 2 ), while the coefficient 2/n is introduced for 
convenience. 

The equations (|3U|l give the optimal coefficients g n : 



.^_l im f ^\(e')\W(e* e ')\Im 

2 r— »1 7 7T 



n > 3. 



W(re i9 ') 

By imposing the constraint (|2*Hj) we find the Lagrange multiplier /i 



(38) 



2(1 - zl) 



c\ff 

go + giZo + 92Z 2 -C(z )f (A K7T )-z* I —\(6')h(6') 

71 



where we used the notation 



h(0) = Im 



,3i0-i* 



1 - z e 



w 



sin[30 - $(0)] - ^ o sin[20 - $ 
1 + Zq — 2z cos 9 



(39) 
(40) 



with 



$(0) = ArgC(e ie ) + 5(6). (41) 

This function is readily obtained using the expression (|2"T|) of the outer func- 
tion C(z) 

HO) = T^ + ^+^Arg[l-e^ + /5(l + e je )]-3Arg[l-e 49 + /? Q (l + e 49 )] +( 5(0) , 

(42) 
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the signs of the first term corespond to 9 > and 6 < 0, respectively. By 
including into the condition J3SJ), the optimal coefficients ([38]) . where /i is 
given in ([39]). and using the Plemelj-Privalov relation 



Ss \ J d "'t^S^t = m + \ J - ^'t^t . (43) 

— 7T — 7T 

where the last integral is a Principal Value, we obtain an integral equation 
for the function X(9) given by 

V{6) ~ X{6) ~ J 2^ X{6 } sin(9-9>)/2 

-(i-^o 2 ) / — mmm, m 

J 71 



where 



2 1 J2 



v{p) = J2 9n sm[ n e-m]-^^ (00 + gizo + g 2 z 2 - c(z )f (A K7T )) h{e). 



n=0 Z 



(45) 

As discussed in Refs. [H], j^Zj, for $(0) given in ([4*2*]) and assuming that the 
physical phase 5(9) is sufficiently smooth (more exactly Lipschitz continuous 
of order greater than 1/2 jEJ), one can show that ([44]) is a Fredholm equation. 

The minimum of the left hand side of ([2*3]) is obtained by inserting the 
solution ([3*5]) into this relation, taking into account the expression (["""*)]) of \i 
and the integral equation (j""""J. We thus obtain the inequality 



I - z 2 f 

0o + 0i + 02 H e - ^ (0o + 0i^o + 02^o - C(z )f (A K7r 

a. 

d9 



+ f -\(6)V(9) <4>"(Q 2 ), (46) 

J 7T 



where \(9) is the solution of the integral equation (JUJ). The inequality ([4*6*]) 
defines the allowed domain for the slope and curvature when both Callan- 
Treiman relation and Watson theorem are imposed. We notice that if the 
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Callan-Treiman relation (JHJ) is removed, the domain allowed for the coeffi- 
cients g , gi and g 2 is given by the inequality [27| 

<?o + 9l + 9l+ I -X(9)V(9) < ip"(Q 2 ) , (47) 

J 71 

where \{6) is the solution of the simpler equation [2*7j 

d9' sin[5/2(0-0') -$(0) + $(0')] 



it-A(^ 

27T 

n=u 



sin(fl - 0')/ 2 

(48) 

Before ending this section we mention that the alternative inequality (fTTj) 
obtained from the dispersion relation (jUJ) leads to similar bounds, which can 
be obtained by replacing in the r.h.s. of the inequalities (J2ZJ), (HHh flSIl) and 
(gZD the quantity ^"(Q 2 ) by (^(Q 2 )/Q 2 )' + ^(0)/Q 4 , and the outer function 
C(z) from (J2U by 

C(z) = OW + + (49) 

where (3q was defined above (|2"Tf . Then, instead of (j2EI), we obtain the 
relations between g n and the Taylor coefficients in (J2J) as: 

<?o = /o(0) 0.00249987, 

gx = / (0) [-0.00154365 + 0.0172731 (r 2 )f w ], (50) 
g 2 = f Q (Q) [-0.00193761 - 0.0452122 (r 2 )f w + 0.00651435 c]. 



3 Results and discussion 

The input required for the numerical evaluation of the bounds is represented 
by: 

1. the function ijj" (Q 2 ) appearing in the right-hand side of the inequalities 

©, flU and ill) 

2. the ratio Fx/F n appearing the Callan-Treiman relation (JHJ), 

3. the phase 6(6) = S^ 2 (t(e l9 )) entering the expression (}4*Tj) of the function 
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Figure 1: Phase shift of the 1=1/2 S-w&ve Ktt amplitude below 1 GeV. 

1 /2 

Solid: phase-shift S (t) calculated in [35J, used in our analysis; dashed: one 

1 /2 PTT^l 

of the curves for 8 (t) given in Fig. 1 of [IH]; dashed-dot: parametrization 
of 5q (0 use d m jSZl- 
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We calculated ip' (Q 2 ) given by the expression (fT3)) taking Q 2 = 4GeV 2 (see 
Ref.|16j). which is sufficiently large to ensure the validity of perturbative 
QCD. At this value of Q 2 the a s correction in ("T3|) is of about 30%, the a 2 
correction is less than 9%, while the higher twist terms contribute with less 
than 1%. The choice Q 2 = 4 GeV 2 is also convenient since predictions for the 
light mi and strange m s quark masses at 2 GeV are now available from lat- 
tice calculations: QCDSF-UKQCD Collaboration [23] predicts m s (2GeV) = 
119(5)(8)MeV and m,(2GeV) = 4.7(2)(3) MeV, where the first quoted er- 
ror is statistical and the second systematic, while MILC Collaboration [2*4] 
predicts m,(2GeV) = 76(0)(3)(7)(0)MeV and m,(2GeV) = 2.8(0)(1)(3)(0) 
MeV, where the errors are from statistics, simulation, perturbation theory 
and electromagnetic effects, respectively. As concerns the strong coupling, 
the value a s (m 2 ) is now well known from hadronic r decays [31], [22]. The 
value given in PDG2004, which we used in our work, is a s {m 2 ) = 0.345 ± 0.03 
[33j, slightly larger than the renormalon-based estimates [51]. We notice that 
the strength of the bounds depend in a monotonous way on the magnitude 
of ip (Q 2 ), a larger value of ip (Q 2 ) at a given Q 2 leads to weaker bounds. 
Therefore, using instead of a s (Q 2 ) at Q 2 = 4 GeV 2 the value a s {m 2 ) is a 
conservative approximation. In this way, we calculate ip (4 GeV 2 ) quite ac- 
curately using the running quark masses and the running coupling at 2 GeV, 
without going through the invariant masses and Aq C d- Introducing stan- 
dard values of the quark condensates [30], we obtained if)" (A GeV 2 ) = 0.00020 
with the quark masses from "23], and ip"(4:GeV 2 ) = 0.000079 with the quark 
masses from |24j . 

The Particle Data Group [33] reports the value 1.22 for the ratio Fx/ 'F n 
appearing the Callan-Treiman relation (J3J. The pseudoscalar decay con- 
stants were recently calculated from the lattice [213, with the result F K /F n = 
1.21 (4) (13), where the first error is statistical and the second systematic. Fol- 
lowing Ref. [""""" . we shall use in our analysis Fk/F^ in the range 1.21 - 1.23. 

1/2 

The last input is provided by the phase shift 5 (t) of the S'-wave Ktt am- 
plitude along the elastic region. We used in our work the phase shift derived 
from a new analysis of Ktt scattering from Roy-Steiner equation [315] and the 
experimental data [3(5]. We took the inelastic threshold at t in = 1 GeV 2 , since 
the inelasticity is equal to 1 up to this energy. In Fig. [TJ we show the phase 
#o it) used in our work, together with the phase shift derived in |18. and 
used in the subsequent analysis [TTj, "TH] of the strangeness changing scalar 
form factors. For completeness we represent also in Fig. ["Ja parametrization 
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Figure 2: Allowed domains for the Kn radius squared and the curvature 
for /o(0) = 0.976 and Fx/F n = 1.21; dashed: unitarity domain given by 
Eq. (|27j) : dashed-dot: domain given by Eq. (|47|) imposing Watson theorem; 
dotted: domain given by Eq. (|32|) imposing Callan-Treiman condition; solid: 
domain given by Eq. ()46|) imposing Callan-Treiman and Watson theorem. 
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of S (t) used in [2IJ for a calculation of the radius (r 2 )^ 71 " based on the 
single-channel Omnes formalism. 

With the input described above we evaluated the inequalities (J27J) . ()32)). 
(I46j) and (|47|) . derived in the previous section. The numerical solution of 
the integral equations (J4"4"|) and (J4"%j) was found easily, since the coefficients 
go, gi and #2 do not appear in the kernel. For a fixed value of /o(0), the 
allowed domain for the slope and curvature is in all cases the interior of 
an ellipse in the plane ((r 2 )f n , c). In order to show the effect produced 
by each additional constraint, we indicate in Fig. |21 the domains given by 
the inequalities (J27|) . (|4*7Jl and (|46|). respectively. The curves shown 

were obtained with /o(0) = 0.976, the quark masses from and Fk/F^ = 
1.21. One can see that the most important effect is produced by the Callan- 
Treiman condition (J2J). Imposing Watson theorem in addition to the Callan- 
Treiman relation leads only to a small reduction of the allowed domain. 

In Fig. El we indicate the allowed domain given by the inequality ()46|). 
which takes into account both Callan-Treiman condition (jSJ) and Watson 
theorem (}3*3*|) . using for the quark masses the values given in Refs. |24 j 
and 23, respectively. As in Fig. El we took / (0) = 0.976 and i^/i^ = 
1.21. Using the quark masses from [24j we obtain the upper bound {r 2 )f n < 
0.26 fm 2 , while with the quark masses from [221 we obtain (r 2 )^ 71 < 0.34 fm 2 . 
In Fig. El we indicate also the prediction of the radius and curvature made 
in by means of an alternative dispersion theory, which gives (r 2 )f n = 
0.192 ± 0.012 fm 2 , c = 0.855 ± 0.051 GeV~ 4 . These values are situated well 
inside the allowed domain derived with both values of the strange quark mass. 
Also, the ;\TT value (r 2 )^ = 0.20 ± 0.05 fm 2 derived in |BJ is consistent with 
the bounds, while the value (r 2 )f n = 0.31 ± 0.06 fm derived in [2ZJ slightly 
violates the upper bound obtained with the s quark mass from [24]. 

Fig. |21 also predicts upper and lower limits on the curvature c for each 
value of (r 2 )f\ For instance, for the ISTRA value j2J: {r 2 )f ±7T = 0.235 fm 2 , 
and the curvature c must be in the range from 0.41 to 0.78 GeV -4 , while for 
the KTeV value |lj : (r 2 )f L7T = 0.165 fm , and the curvature must be larger 
than 0.67 and smaller than 1.23 GeV -4 . Thus, the second order term in 
the Taylor expansion (0) seems to be important and can not be neglected in 
experimental fits. The allowed range of c given in Fig. |2]for each value of the 
slope may be useful as an additional constraint in the representation of the 
experimental data. Up to now we indicated the allowed domains obtained 
with the form (j2J) of the dispersion relation. The relation © leads to similar 
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Figure 3: Allowed domain for the Kn radius squared and curvature, given by 
Eq. ()46)1 which incorporates Callan-Treiman relation and Watson theorem. 
Dashed: values obtained with m s = 119 MeV, m u = 3 MeV cf. |2H|; solid: 
values obtained with m s = 76 MeV, m u = 2.8 MeV cf. [21]. The diamond is 
the determination made in (TH] for the same input values /o(0) = 0.976 and 
Fk/F-k = 1-21, while the circle with error bars is the final determination in 
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<r > s (fm ) 

Figure 4: Allowed domain for the K7r radius squared and curvature given by 
Eq.f^HI) using m s = 76 MeV, m u = 2.8 MeV: solid: domain obtained with the 
dispersion relation dotted: domain obtained with the dispersion relation 
©. The points are the determination made in [19] . 
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results. In this case, in the r.h.s. of the inequalities ()27j). ()32jl. (|47)) and (|35j) 

appears the quantity (ip{Q 2 ) /Q 2 )' + ^(0)/Q 4 . We evaluated this quantity at 
Q 2 = 4GeV 2 using the QCD expression (JT3)l and the inequality ((7|), which 
leads to the values 0.00022 with the quark masses from [21], and 0.00038 
with the quark masses from [23]. In Fig. 01 we present for comparison the 
allowed domain for the radius and curvature obtained with the two forms of 
the dispersion relation discussed above. 

It is of interest to combine the bounds on the Taylor coefficients of fo(t) 
derived in the present work with the predictions of x?T. The expression of 
the scalar form factor at order p 6 in x?T is jH] 

f (t) = F + (0) + A(t) + Fk/ ^ ~ 1 1+ JL(2C r 12 + C r u )Z K7r t- ±Cy*, (51) 

where 

f + (o) = 1-0.008- ^(q 2 + q 4 )AL, 

A{t) = -0.259t + 0.84t 2 + 1.291t 3 , (52) 

where the chiral constants are taken at the scale \x = M p [S], PH]- The 
inequalities derived in the previous section can be expressed as an allowed 
domain for the chiral constants C[ 2 and C r M . As in Ref.JTH] we use as inde- 
pendent parameters C[ 2 and the sum C[ 2 + C£ 4 . For Q 2 = 4GeV 2 , the first 
three Taylor coefficients g n defined in (124)1 are expressed in terms of these 
parameters as 

= 0.00163 - 9.23 (C r 12 + C r M ), 
gx = _o.000104 + 77.2 C[ 2 + 87.84 (C[ 2 + C r u ) . (53) 
g 2 = -0.00226 - 713.3 C[ 2 - 241.93 (C r 12 + C r M ). 

Then inequality ()46)) gives the allowed domain shown in FigEl where we 
indicate also the result derived in ^2] and the recent determination of the 
chiral low-energy constants in [3*H| . 

From Figs. EE it is seen that the predictions made in are much more 
precise than the unitarity bounds derived by us, even with the implementa- 
tion of Callan-Treiman and Watson theorem. It is of interest to see what is 
the additional input used in JH] , responsible for the very small uncertainties 
quoted for their results. In Ref. ^7] the authors solve the coupled channel 
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Figure 5: Allowed domains for the low energy constants C{ 2 and C{ 2 + 
using Eq. (|46|). with the quark masses from [21] and Fk/F^ = 1.21. The 
circle with errors is the determination in JH] , the square is the prediction in 
[38J, where the quoted errors represent the variation with the scale /i in the 
range from M, q to 1 GeV. 
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unitarity equations for the Kn, Krj' and Kt] form factors, using the chiral 
resonance model and various JC-matrix parametrizations of the scattering 
amplitudes. By construction, the Watson theorem is satisfied by the Kn 
form factor fo(t), while the two integration constants entering the solution 
of the coupled channel equations are fixed by the values of fo(t) at t — and 
t = A Kw cf. Callan-Treiman relation (jHJ). Therefore, f (t) calculated in [T7] 
satisfies the same low energy constraints which we imposed. However, in our 
formalism the form factor is constrained above the inelastic threshold only by 
the L 2 condition (or (fTTj)). This puts a weak restriction on the modulus 
|.F(t)|, which is allowed to increase like t at asymptotic energies, while in 
[T7j the 1/t asymptotic decrease predicted by perturbative QCD is imposed. 
Analyticity is also implemented in a different way in the two approaches: in 
our formalism it is expressed by the relations (J23j) and ([24)1. where the real 
coefficients g n obey only the constraints that were explicitly incorporated. 
So, we work with the most general class of functions which satifies these 
constraints. On the other hand, in ^7] the form factor is expressed by a 
dispersion relation with the spectral function related to the meson - meson 
transition amplitudes, described by several specific parametrizations. Choos- 
ing a restricted class of functions leads to a smaller error, but may introduce 
a bias, so the errors given in J7| might be underestimated. 

Of course, the information provided by the inelastic channels is impor- 
tant. The contribution of other two-particle states can be implemented in a 
rigorous way in the present formalism, as shown in [27| (see also |2H|). The 
higher states in the unitarity relation (JSJ) involve the modulus squared of 
the Kt] and Krf form factors above the corresponding unitarity thresholds. 
Their unphysical cuts can be eliminated by exploiting Watson theorem with 
a generalized Omnes formalism j2Zj. We shall investigate this problem in a 
future work. Other generalization, easily implementable, is the introduction 
of higher derivatives at zero momentum transfer, which are related to the 
chiral constants of %PT. 

We mention finally that the inequalities derived above can be used also 
in the opposite way, for bounding the quark masses using the information 
on the low energy expansion parameters of the strangeness changing scalar 
form factors. This problem was investigated for instance in jlU] and |4"T] . 
Our approach offers the possibility to include simultaneously various pieces 
of information, like Callan-Treiman constraint Q and Kir phase shift 5^ 2 
through Watson theorem. For instance, using the predictions of the radius 
and the curvature made in ^1 we obtain the lower bound m s (2 GeV) — 
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m u (2 GeV) > 38.6 MeV. The prediction made in jH] based on a QCD sum 
rule is tighter since the authors use the spectral function of the two point 
function ip calculated from the values of the Rn, Kr\ and Krf form factors 
on the unitarity cut. 

4 Conclusions 

We have derived and evaluated unitarity bounds on the slope and curvature 
of the scalar strangeness changing form factor fo(t). The method incorpo- 
rates in an exact mathematical way low energy theorems like Callan-Treiman 
condition and Watson final state interaction theorem. Our results indicate 
that the curvature can not be neglected in the representation of the exper- 
imental data on decay. Using the calculation of the form factor to two 
loops in %PT, we expressed our results as an allowed domain for the chiral 
constants C{ 2 and C[ 2 +C^. Our bounds confirm the values obtained recently 
in an coupled channel dispersion approach, but are much weaker (however 
the uncertainties given in [F3\ might be underestimated). Generalizations of 
the method, including the Kt] and Krf form factors in the unitarity sum for 
the spectral function, and their higher derivatives at the origin, can increase 
the predictive power of the present formalism. 
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